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ABSTRACT 

Broad infrared emission featnres (e.g., at 3.3, 6.2, 7.7, 8.6, and 11.3 /im) 
from the gas phase interstellar medium have long been attributed to polycyclic 
aromatic hydrocarbons (PAHs). A signihcant portion (10 to 20%) of the 
Milky Way’s carbon reservoir is locked in PAH molecules, which makes their 
characterization integral to our understanding of astrochemistry. In molecular 
clouds and the dense envelopes and disks of young stellar objects (YSOs), PAHs 
are expected to be frozen in the icy mantles of dust grains where they should 
reveal themselves through infrared absorption. To facilitate the search for 
frozen interstellar PAHs, laboratory experiments were conducted to determine 
the positions and strengths of the bands of pyrene mixed with H 2 O and D 2 O 
ices. The D 2 O mixtures are used to measure pyrene bands that are masked by 
the strong bands of H 2 O, leading to the hrst laboratory determination of the band 
strength for the CH stretching mode of pyrene in water ice near 3.25 pm. Our 
infrared band strengths were normalized to experimentally determined ultraviolet 
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(UV) band strengths, and we find that they are generally ~ 50% larger than 
those reported by Bouwman et ah (2011) based on theoretical strengths. These 
improved band strengths were used to reexamine YSO spectra published by 
Boogert et al. (2008) to estimate the contribution of frozen PAHs to absorption 
in the 5 to 8 /rm spectral region, taking into account the strength of the 3.25 /rm 
CH stretching mode. It is found that frozen neutral PAHs contain 5 to 9% of the 
cosmic carbon budget, and account for 2 to 9% of the unidentified absorption in 
the 5 to 8 /im region. 


1. Introduction 

Unidentified infrared bands discovered almost 40 years ago in the interstellar medium 
(ISM) are now attributed to the CH and CC vibrational modes of polycyclic aromatic 
hydrocarbons (PAHs). They are estimated to account for 10 to 20% of the total carbon 
reservoir in the ISM (Peeters 2011). Absorption features at 3.25 (Sellgren et al. 1994; Brooke 
et al. 1996, 1999), 6.2 (Keane et al. 2001), and 11.3 /im (Bregman et al. 2000) have been 
attributed to PAHs in a handful of YSO spectra, while PAH emission bands are weak or 
absent toward embedded YSOs (Geers et al. 2009). In molecular clouds and in the disks 
and envelopes of embedded young stellar objects (YSOs), PAHs, along with other gas phase 
species, are expected to freeze out on the icy mantles of dust grains. In the dense ISM, 
ice mantles are important sites for prebiotic chemistry, where the long residence times of 
condensed molecules give them the best opportunity for interaction. Indeed, laboratory 
experiments have shown the ease of producing biomolecules via irradiation and heating of 
interstellar ice analogs (e.g. Bernstein et al. 2002; Munoz Caro et al. 2002; Kaiser et al. 2013). 
Much work has been done to understand PAH chemistry in ice matrices (e.g. Gudipati & 
Allamandola 2003; Gudipati 2004; Gudipati & Allamandola 2004). In particular, PAHs 
frozen in amorphous H 2 O ice ionize quickly when irradiated by UV light. Larger PAHs 
remain stable up to 120 K, at which point the ice crystallizes and PAHs begin to react with 
the ice matrix itself (Gudipati & Allamandola 2006; Allamandola 2011), while smaller PAH 
ions react at lower temperatures to form hydrogenated and oxygenated complex molecules 
(Gudipati & Yang 2012). 

Significant attention has also been given towards the task of quantifying PAHs in the 
ISM. An extensive online database at www.astrochem.org (Bauschlicher et al. 2010; Boersma 
et al. 2014) contains both theoretical and experimental spectra of PAHs matrix-isolated in 
argon. These spectra greatly enhance our understanding of PAH signatures, but spectra of 
PAHs embedded in interstellar ice analogs are necessary to truly constrain PAH abundances 
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within the dense ISM. Experiments have been carried out to determine band strengths for 
both neutral (Bernstein et ah 2005; Bouwman et al. 2011) and ionized (Bernstein et ah 2007) 
pyrene frozen in water ice, but the absolute infrared (IR) band strengths reported in these 
cases were normalized to theoretically calculated band strengths. Sometimes, the region of 
interest reported was constrained between 1650 to 1000 cm“^ (6 to 10 pm) because the H 2 O 
ice absorption there is less signihcant and more linear than in other regions, and none of 
these studies address the CH stretching mode at 3.25 pm, which is most easily accessible 
with ground-based telescopes. Band strengths for naphthalene frozen in water ice have been 
reported as well (Sandford et al. 2004). 

The aim of our work is to improve upon previous band strength estimates by recording 
PAH spectra in both infrared and ultraviolet-visible (UV-Vis) bands, then normalizing the 
results to previous, direct measurements of UV absolute band strengths (Berlman 1971; 
Dixon et al. 2005). This paper reports our results for pyrene (CigHio) frozen in either H 2 O 
or D 2 O ice. Pyrene was chosen as a representative PAH molecule because it has been widely 
studied in the laboratory and a signihcant amount of laboratory data exists. The pyrene 
radical cation is proposed to be the carrier of some of the diffuse interstellar bands (DIBs) 
(Salama & Allamandola 1992). Additionally, pyrene is the smallest available PAH that is 
convenient to handle but at the same time compact with a C/H ratio of 1.9 which is closer 
to the larger PAH C/H ratio expected in the interstellar medium (see Section 4.3). We 
included D 2 O in our experiments because its absorption features are redshifted compared 
to those in H 2 O, allowing easier measurement and detection of some pyrene bands which 
would normally be drowned out by H 2 O features. With revised band strengths in hand, we 
attempted to identify the previously reported 3.25 pm band (Sellgren et al. 1994; Brooke 
et al. 1996, 1999) in the set of YSO spectra published by Boogert et al. (2008) and, where 
found, constrained PAH column densities and the contribution of PAH absorption to the 5 
to 8 pm absorption region. Section 2 describes the experimental setup for obtaining spectral 
measurements while Section 3 details our procedure for calculating absolute band strengths. 
Section 4 describes the astrophysical implications for quantifying PAHs in the dense ISM. 
Section 5 summarizes our findings and suggests possible trajectories for further expansion of 
this work. 


2. Experimental Setup 

The setup for our laboratory experiment was very similar to that described in Barnett 
et al. (2012); a schematic is shown in Figure 1. In summary, the chamber for this experiment 
was cooled to 30 K at which temperature the base pressure in the chamber was in the high 
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10“® to low 10“® mbar. Water vapor was passed over pyrene crystals which were heated to 
~ 55 °C for sublimation under vacuum. This mixture of sublimated pyrene and water vapor 
was directed into the chamber, where it subsequently froze onto a potassium bromide (KBr) 
window. During the ice deposition, pressure in the chamber was typically around 10“^ to 
10“® mbar due to the water vapor lead into the chamber. Under these conditions, the pyrene 
concentration was less than three percent in all samples, ensuring that pyrene molecules 
were adequately isolated from each other within the ice matrix. The pyrene concentration 
was further monitored and controlled as necessary during the experiments by measuring the 
ice him thickness, adjusting the pyrene sublimation temperature, and monitoring the IR 
and UV-Vis absorption features. Ice film deposition rates were approximately 1 pm 
per honr. Transmission spectra throngh samples ranging in thickness from a 
few hnndred nanometers to 5 pm were collected via Fourier transform infrared (FTIR) 
spectroscopy and a UV-Vis spectrograph using a deuterium/halogen lamp as a continuum 
source. The spectral resolution and spectral coverage of the FTIR spectrometer were 1 cm“^ 
and 8000 to 400 cm“^, respectively. The UV-Vis spectrometer had a resolution of ~ 0.75 nm 
in the 210 to 1100 nm wavelength range. Each instrument was connected to a separate 
port of the octagonal chamber, allowing IR and UV spectra to be collected consecutively. 
The optical window was oriented in such a way that it was as close to perpendicular as 
possible to the deposition side, at the same time allowing both FTIR and UV-Vis spectral 
measurements without rotating the cryostat on which the optical window was mounted. This 
conhguration gave us highly reproducible spectra as we grew the ice slowly on the optical 
window. Roughly half the experiments were conducted with H 2 O and half with D 2 O in 
order to recover better measurements in areas with significant ice absorption. During D 2 O 
matrix experiments, even though we passivated the chamber with D 2 O, there was still, at 
times, obvious contamination of H 2 O ice features in the spectra, necessitating great care 
when interpreting band strengths in regions of overlap between ice features. 

Atmospheric water vapor in the path of the FTIR instrument and the detector was 
minimized by passing gaseous nitrogen through the enclosed regions. However, there were 
several openings, which were necessary to optimize the instrumentation, that could not be 
completely sealed off. As a consequence, vapor-phase water absorption is always present in 
our infrared spectra. Fortunately, due to their characteristic absorption, we have routinely 
subtracted the water vapor absorption from the ice spectra, as will be discussed in the next 
section. Water vapor spectra were derived from measuring the spectra of the experimental 
set-up at room temperature before and after passing dry nitrogen gas at a very high flow rate 
(10 times more than the normal flow rate of ~ 20 L/min), eliminating a signihcant amount 
of water vapor over a short period of time. 
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Deuterium/Halogen Lamp 
(UV-Vis Light Source) 



Fig. 1. — Experimental setup. The pressure in the chamber was kept at ~ 10“^ mbar, 
and the temperature of the KBr window was kept at 30 K. The frozen pyrene sample was 
generated by passing water vapor over sublimating pyrene crystals. This gaseous mixture 
then froze onto the KBr window. Transmission spectra were recorded using FTIR and UV- 
Vis spectrometers. 


3. Analysis 

Absolute IR band strengths of ices are difficult to calculate directly from laboratory 
experiments, because column densities require a determination of the sample density, thick- 
ness, and mixing ratio. In order to determine absolute IR band strengths of pyrene frozen 
in water ice, we correlated our measured IR pyrene bands to the hrst electronic transition 
of pyrene. The absolute band strength of this UV band is easily calculated from previously 
published data. 
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3.1. UV Band Strength And Integrated Absorption 


The molar absorptivity, e (M“^ cm“^), of the hrst electronic transition of pyrene (in 
cyclohexane) was measured by Berlman (1971). This transition occurs in the UV frequency 
range of ~ 35000 to 28800 cm“^ (286 to 347 nm). Dixon et ah (2005) duplicated the 
experiment and published their results in the online database PhotochemCAD^. Their results 
were normalized at the peak of the hrst electronic transition to the value e = 54000 cm“^ 
listed by Berlman (1971), which agrees with the results found by Thony & Rossi (1997) for 
gas phase pyrene at 150 K. From these measurements, we directly calculated the inte- 
grated band strength of the hrst electronic transition of pyrene using Equations (l)-(4). 
Althongh these experiments measured pyrene in cyclohexane rather than in ice, 
the total oscillator strength of any given electronic transition remains constant 
(a molecular property), irrespective of the solvents used, unless strong elec- 
tronic interactions between the solvent and the PAH molecule occurs that could 
change electronic and vibrational transitions significantly. This is not the case 
in cyclohexane which is a non-polar solvent. Amorphous ices behave like non- 
polar solvents with electronic polarizability similar to argon matrix, so strong 
interactions between pyrene and amorphous water ice are negligible (Gudipati 
2004). Therefore, the assumption that the integrated UV band strength for the 
first electronic transition of pyrene is equivalent in amorphous water ice and 
cyclohexane is valid. 

Integrated band strengths. A, are calculated using 


r(x)dx 2.303 a{x)dx 

“ ii 


( 1 ) 


where r = ln{Io/I) is the optical depth of the feature, x is the wavenumber (1/A, cm“^), 
N is the column density (cm“^) of the material producing the feature, and a = logio{Io/I) 
stands for absorbance, which is the preferred unit of experimental scientists. 

After some manipulation, the Beer-Lambert law (2) can be substituted into (1) to 
determine the band strength directly from molar absorptivity. In Equations (2)-(4), c is 
molar concentration (1 M = 0.001 mol cm“^), I is path length (cm), n is number density 
(cm“^), and Na is Avogadro’s constant (6.02 x 10^^ mol~^). 


a = eel 


( 2 ) 


^ WWW. photochemcad.com 
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N = nl 


(^ 0.001 


mol cm 

w~ 


NacI 


( 3 ) 


A 


2.303 e(x)dx 


0.001 


mol cm 
M 


Na 


( 4 ) 


Before performing any of the following analysis, raw spectra were smoothed using an 
FFT filter (Bowman 2006) with a window^ of p = 5 in order not to compromise actual 
signal. The integration of molar absorptivity was carried out after subtracting a local, 
linear baseline anchored at 35000 and 28800 cm“^. We calculated the UV band strength as 
2.578 X 10“^^ cm/molecule. 

The errors in molar extinction coefficients come not from the spectra, but determining 
the weight of the substance in solution (molarity) as well as whether the molecules form 
aggregates in solution or stay as monomers (Sin & Duhamel 2008). A good comparative study 
on absorption of pyrene in various solvents is given in Ray et ah (2006), where significant 
changes in peak positions are seen. However, their integrated molar absorptivity remains 
essentially unaltered, because if the height decreases, the bandwidth increases. Hence, we 
estimate that the experimental band strengths derived for pyrene UV spectra are very 
accurate within a few percent uncertainty at the maximum. 

Integrated absorbances of the first electronic transition of pyrene were calculated from 
our spectra using the same linear baseline as described above. Subtracting a local baseline 
is important for frozen pyrene, because absorption due to the ice matrix itself is significant. 
Errors on the integrated UV absorbances were determined by shifting the baseline by 3 nm 
in either direction. In the rest of this work, the pyrene column densities calculated from the 
UV spectra were used to derive band strengths for the infrared transitions. 


^See http://wiki.originlab.coin/~originla/wiki2/index.php?title=X-Function:Smooth. Note, however, a 
typo in the equation listed there. The term (1 — in the last line of the equation should instead be {n — 


Wavelength (lim) 

0.278 0.286 0.294 0.303 0.3T2 0.323 0.333 0.345 0.357 



36000 35000 34000 33000 32000 31000 30000 29000 28000 

Wavenumber (cm“^) 


Fig. 2. — Panels A and B show intermediate steps for obtaining the absolute UV band 
strength for the hrst electronic transition of pyrene. An absorption spectrum selected 
from our sample of pyrene in H 2 O is shown in blue. The molar absorptivity spectrum 
of pyrene in cyclohexane, published by Dixon et ah (2005), is shown in red. Panel A shows 
original spectra, without baselining. Panel B shows the same spectra, after subtracting linear 
baselines anchored at 35000 and 28800 cm“^ (dotted lines). 


3.2. IR Integrated Absorption 

To determine the strength, peak position, and width of the IR pyrene features, contam- 
ination by H 2 O vapor and baseline curvatures (due to ice absorption) need to be removed 
hrst. Contamination by water vapor absorption is signihcant in many of the spectra from 
~ 2025 to 1275 cm“^ (see Figure 3). Without removal, such contamination greatly affects 
the identihcation and analysis of almost half of the IR pyrene bands. To obtain a water 
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vapor template, we took the ratio of two background spectra. The water vapor residual 
was cleaned by the same FFT smoothing described above, and a local, low-order polynomial 
baseline was subtracted. The amplitudes of the hve largest water vapor peaks unaffected 
by ice absorption were measured in each spectrum, and the water vapor template, scaled to 
those amplitudes, was then subtracted to remove the contamination. 


Wavelength (tim) 

6.061 6.452 6.897 7.407 



Wavenumber (cm 


Fig. 3. — Typical example of water vapor subtraction. Regions where pyrene features exist are 
marked in grey. The original spectrum is shown in black, and the corrected spectrum (green) after 
subtraction of water vapor contamination (blue) is offset below. 


To determine the integrated absorption of pyrene features in the IR, a global baseline 
to remove ice absorption was determined through implementation of the following method. 
First, pyrene features were masked (see Table 1), and the remaining spectrum was smoothed 
using Gaussian convolution. Then, the smoothed spectrum was interpolated across feature 
regions using a second-order polynomial. This global baseline was subtracted from the 



10 


original spectrum to obtain a residual pyrene spectrum for which position, FWHM, and 
integrated area was determined for each band. Figure 4 provides an example of this method. 

The error in the integrated area of IR bands was calculated by taking into account the 
standard deviation of absorption in the pyrene residuals. Since the noise can vary slightly 
with wavenumber, the standard deviation in absorption for a given feature was calculated 
within a region ± 3 times the width of the masking region for that feature, ignoring any 
other pyrene features within these limits. 



3100 3000 

Wavenumber {cm"^) 


Wavelength (iim) 

6.061 6.452 6.897 7.407 8.000 8.696 9.524 



Wavenumber (cm 



950 850 750 650 

Wavenumber (cm'^) 


Fig. 4. — Typical examples of pyrene/H20 and pyrene/D20 spectra, after water vapor removal, 
are plotted with their baselines in the regions of pyrene absorption. H 2 O is in blue, and D 2 O 
is in green. The PAH residual spectra are offset below. Due to the way in which features were 
masked during baselining, not all PAH features were measured in each spectrum. Features that 
were measured are marked with a filled circle, and features that were ignored are marked with a 
cross. Note that the central panel has a much smaller y-axis scaling. 


3.3. Absolute IR Band Strengths 

The absolute IR band strengths were determined by correlating a strong infrared band 
with the calibrated UV band. This was challenging because the UV bands of pyrene tend to 
saturate fairly quickly while the weaker IR bands can take much longer to build up enough 
signal-to-noise for measurement. Figure 5 shows the correlation of the 1436 cm“^ band, 
as shown in Figure 4, with the UV band. This band is both fairly strong and relatively 


Table 1. Masking Regions for Global Continuum 


Band Identihcation 

Masking Region (cm ^) 

A 

3078 to 

3032 

B 

1607 to 

1584 

C 

1493 to 

1483 

D 

1477 to 

1463 

E 

1460 to 

1444 

F 

1441 to 

1426 

G 

1413 to 

1402 

H 

1321 to 

1308 

I 

1251 to 

1235 

J 

1194 to 

1171 

K 

1103 to 

1091 

L 

1072 to 

1059 

M 

864 to 

837 

N 

826 to 

817 

0 

767 to 

746 

P 

722 to 

700 


Note. — For this analysis, features were 
manually selected by visual inspection of the 
laboratory spectra using Bouwman et ah (2011) 
and the experimental results for pyrene in argon 
from www.astrochem.org (Bauschlicher et ah 
2010; Boersma et al. 2014) as a guide. Although 
an additional feature at 1136 cm“^ was reported 
by Bouwman et al. (2011), it was not detected 
above 3a in our sample. 
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unaffected by ice absorption in either H 2 O or D 2 O ice. Since IR and UV spectra were 
taken consecutively rather than concurrently, IR measurements had to be interpolated to 
the observation times of the UV spectra. After determining the absolute band strength 
for the 1436 cm“^ band, the strengths of all bands were determined by correlating them 
with each other. Linear correlations were determined using the IDL program MPFITEXY 
(Williams et ah 2010) which utilizes the MPFIT package (Markwardt 2009). MPFIT is a 
curve htting program that uses a robust, non-linear least squares method. MPFITEXY was 
used to take into account x and y error bars and force the y-intercepts through the origin. 
Table 2 lists the position and FWHM of each band averaged over all spectra along with the 
hnal calculated band strengths. 



Fig. 5. — Correlation between the 1436 cm“^ IR band and the 35000 to 28800 cm“^ UV 
band. Triangles show H 2 O samples, and squares show D 2 O samples. Small circles show how 
the correlation breaks down as the UV band begins to saturate. Non-linear growth of the UV 
features becomes signihcant above a peak absorbance of 0.8 (or r of 1.84). Only unsaturated 
data points were used to determine the best ht linear correlation with the condition that it 
passes through the origin. 


Table 2. IR Pyrene Band Characteristics 
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3.3.1. Comparison of Pyrene Band Strengths in H 2 O and D 2 O Ices 

Spectra were recorded for pyrene in D 2 O as well as in H 2 O ice in order to recover bands 
masked in laboratory spectra by the strong H 2 O features (e.g., the 3.25 /rm CH stretch 
mode). This requires an assumption that CH stretching modes of pyrene behave identically 
whether the pyrene is frozen in H 2 O or D 2 O ice. Figure 6 demonstrates that this assumption 
is valid at least for the CC stretching and CH in plane bending modes of pyrene. It is not 
clear that this is the case for the CH out of plane bending modes. Unfortunately, 
beyond 10 /rm, ice absorption in both H 2 O and D 2 O is quite strong and makes 
characterization of pyrene bands in that region more difficult. Local baselines 
could not be adequately determined for many of the pyrene bands beyond 10 /rm, 
but the 850 cm ^ band appears to be measurable in both ice matrices. This was 
the only band in our analysis which was measured in both ices and gave different 
results depending on the ice matrix, as noted in Table 2. 


Wavelength (nm) 

5.882 6.250 6.667 7.143 7:692 8.333 9.091 10.000 



Wavenumber (cm ') 


Fig. 6. — Comparison of the pyrene residual spectrum averaged across H 2 O (blue) and D 2 O 
(green) experiments. The average spectra were normalized at 1436 cm“^. Regions of obvious 
pyrene features are marked below in red. Some regions had to be excluded from the analysis 
due to inability to accurately determine baseline ice absorption, e.g., near 1200 cm“^ in D 2 O 
ice. 


3.4. Comparison With Previous Laboratory Work 

Figure 7 compares our results with those published by Bouwman et al. (2011). While 
the positions of features measured in both works are essentially identical, the band strengths 
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we report are ~ 50% larger than those reported by Bouwman et ah (2011). Apparently, 
the theoretical band strength calculations used by Bouwman et ah (2011) consistently 
underestimate the values measured by our UV /IR correlation method. 




Fig. 7. — Peak position (left panel) and band strength (right panel) comparison between 
this work and that of Bouwman et ah (2011). The error bars along the y-axis in the left 
panel are smaller than the plot points. Our calculated band strengths are roughly 1.5 times 
those published by Bouwman et ah (2011). 


4. Astrophysical Implications: PAH Absorption Features in YSO Spectra and 

the Carbon Budget 

As illustrated here for pyrene, PAH IR spectra are characterized by a CH stretching 
feature near 3.25 microns, CC stretching and CH in plane bending modes between 5 and 
10 microns, and CH out of plane bending modes between roughly 10 and 15 microns. The 
detection of a 3.25 pm band, attributed to the CH stretching mode of PAHs, was reported 
for a handful of YSO spectra (Sellgren et ah 1994; Brooke et ah 1996, 1999). We looked 
for the presence of such a feature in the YSO spectra published by Boogert et ah (2008). 
For those spectra where a feature was detected, we attempted to quantify the PAH column 
density and the contribution of absorption by PAHs to the 5 to 8 pm region. Since we only 
have laboratory data for pyrene, the following measurements are reported as a function of 
the number of CH and CC bonds present. This of course requires the assumption that band 
strengths per bond remain relatively constant amongst different PAH species. Figure 8 shows 
the validity of this assumption. Here, experimentally determined band strengths per bond for 
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neutral PAHs in argon from the www.astrochem.org database are shown as a function of the 
number of CH and CC bonds per molecule. For the CH stretching and CC stretching/CH 
in plane bending modes, band strengths were integrated from 3200 to 2800 cm“^ and 2000 
to 1250 cm“^, respectively. 




10 20 30 40 50 60 

Number of CC Bonds 


70 


Fig. 8. — Experimentally measured band strengths for the CH stretching mode near 3.25 /rm 
(left) and the CC stretching/CH in plane bending modes between 5 to 8 /rm (right) from the 
www.astrochem.org database (Bauschlicher et ah 2010; Boersma et ah 2014) for neutral PAH 
species in argon (circles). Band strengths per bond remain fairly constant across reported 
PAH species, although the 5 to 8 /rm band strengths do change mildly with a slope of 
0.09 km mol~^ per CC bond per total number of CC bonds (excluding the outlier due to 
C 4 sH 2 o)- Dashed lines show the best htting correlations between band strengths per bond 
and number of bonds, while solid lines show the best htting constant band strengths per 
bond. Stars mark the band strengths per bond for pyrene in water ice. 


4.1. 3.25 yum Feature 

We analyzed the 3.25 /im absorption band in the YSO spectra published by Boogert et ah 
(2008). In summary, a total of 21 objects were observed with sufhcient wavelength coverage 
to examine their spectra for a feature near 3.25 /rm. Six of those objects were massive YSOs 
observed with the Short Wavelength Spectrometer on board the Infrared Space Observatory 
(also published by Gibb et al. 2004). The others were low mass YSOs observed in the L-band 
with ground-based instruments and longward of 5 fim with the Spitzer Space Telescope. 

Ground-based spectra were much noisier than the ISO spectra, so they were hrst 
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smoothed with the boxcar technique using a window of 5 points. Only a subset of the ISO 
spectra analyzed included flux density error bars, which were subsequently propagated to 
determine error bars for the optical depth spectra. To determine error bars for the remaining 
optical depth spectra, a line was £t to a region ±0.025 /rm on either side of each data point, 
and the standard deviation of that data set with respect to the linear ht was assigned as the 
error bar value for the relevant point. 

For each spectrum, a low-order (3'^^ to 5**^) polynomial was £t to a continuum region, 
dehned in Table 3, depending on the absorption features of each spectrum. Next, the residual 
spectra were fit with up to five Gaussians, located at ~ 3.25, 3.32, 3.36, 3.47, and 3.53 /rm, 
using MPFIT and MPFITFUN (Markwardt 2009). The 3.47 /im feature, often seen in YSO 
spectra, has been attributed to solo hydrogens attached to “diamond” -like carbon clusters 
(Allamandola et ah 1992) or to ammonia hydrate (Dartois & d’Hendecourt 2001). The 
remaining model features at 3.32, 3.36, and 3.53 /im correspond to the expected locations 
of solid CH4, C 2 H 6 , and CH3OH, respectively (e.g., Whittet 2003). In the ground-based 
spectra, the 3.32 /rm feature is heavily contaminated by telluric CH4 lines. 

Of the ISO spectra, only three objects have a significant (> 3cr) feature near 3.25 /rm. 
Brooke et ah (1999), using ground-based spectra, listed previously reported measurements of 
the feature for Mon R2 IRS 3 and S140 IRS 1, but they excluded measurement of an apparent 
feature at 3.2 /im for GL 2136 as they claim that its coincidence with the inflection point 
of crystalline ice absorption interferes with accurate measurement. These objects are shown 
in Figure 9, and the parameters of their model fits are listed in Table 4. The integrated 
optical depth measurement of the 3.25 /rm feature for Mon R2 IRS 3 is in good agreement 
with previous results by Sellgren et ah (1994, 1995), and our measurement for S140 IRS 1 is 
similar to that reported by Brooke et ah (1996) to within a factor of two. 

For our ground-based YSO spectra, in order to increase the signal-to-noise ratio, all 
spectra were averaged to recover a global 3.25 /rm feature. The average, continuum sub- 
tracted spectrum, binned to 0.002 /rm, is shown in Figure 10 on an optical depth scale. The 
best £t Gaussian to the region of interest is located at 3.258 ± 0.004 /im with a maximum 
optical depth of 0.021 ± 0.006 and a FWHM of 24.13 ± 6.93 cm-^ (0.026 ± 0.007 /rm). 


Table 3. Continuum Definition 


Object 


R.A. & Dec. 
(J2000) 

Continuum Regions 
(^m) 

Mon R2 IRS 3 

06’"07“ 

A7".80 -06°22'55".0 

3.15 to 3.20 

3.28 to 3.30 

3.70 to 3.80 

S140 IRS 1 

22^19“ 

T8M7 -h63°18'47".6 

3.15 to 3.20 

3.28 to 3.36 

3.70 to 3.80 

GL 2136 

18^22“ 

26".32 -13°30'08".2 

3.14 to 3.17 

3.25 to 3.33 

3.70 to 3.80 

ground-based spectra 



3.15 to 3.20 

3.28 to 3.36 

3.70 to 3.80 
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MonR2 IRS3 



S140 IRS1 



GL2136 



Fig. 9. — Optical depth spectral models after continuum subtraction. The model 
fits, composed of up to 5 Gaussians, are shown with solid blue lines. Each 
component of the model is shown separately with dashed red lines. 
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Fig. 10. — Average residual spectrum from ground-based observations, binned to 0.002 pm 
(red). The individual residual spectra after continuum subtraction are all overplotted in 
grey. A Gaussian, fit to the prospective PAH feature, is shown in blue. The sharp feature 
near 3.32 pm is due to telluric CH4. 


4.2. 5 to 8 pm Region 

Boogert et ah (2008) found that the 5 to 8 pm absorption regions in YSO spectra are 
dne to at least hve distinct components in addition to the bending mode of H 2 O ice, each 
of which likely have multiple carriers. Since PAHs exhibit a number of features in this same 
region due to CC stretching and CH in plane bending modes, we hypothesized that some of 

Table 4. 3.25 pm Absorption Feature 


Object 

A 

FWHM 

'^max 


pm 

fim 

cm ^ 


Mon R2 IRS 3 

3.244 ± 0.002 

0.052 ± 0.004 

49.3 ± 3.4 

0.046 ± 0.003 

S140 IRS 1 

3.247 ± 0.001 

0.046 ± 0.003 

44.1 ± 3.1 

0.028 ± 0.002 

GL 2136 

3.206 ± 0.003 

0.051 ± 0.006 

49.4 ± 6.0 

0.163 ± 0.018 

ground-based average 

3.258 ± 0.004 

0.026 ± 0.007 

24.1 ± 6.9 

0.021 ± 0.006 
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the absorption could be due to PAHs frozen in ice mantles. 

Assuming that any 3.25 /mi feature present is due to neutral PAHs frozen in ice, we used 
our previous calculations of absolute band strengths for pyrene to constrain the potential 
PAH contribution to the 5 to 8 /tm absorption region (fpAH, see Table 5), after subtraction of 
the H 2 O bending mode. Ours is only a first order approximation since the PAHs present in 
the ISM are likely to consist of a mixture of many species, as clearly illustrated by Figure 11. 


Wavelength (am) 
3.125 3.279 3.448 



Wavenumber (cm ') 


Wavelength (am) 

4.762 5.128 5.556 6.061 6.667 7.407 8.333 9.524 11.111 13.333 16.667 



2100 1950 1800 1650 1500 1350 1200 1050 900 750 600 

Wavenumber (cm"') 


Fig. 11. — The observed absorption spectrum (grey) of Mon R2 IRS 3 is compared to the 
laboratory spectrum of pyrene. The left panel shows the residual 3.25 /im feature after 
subtraction of a local baseline (see Section 4.1). The YSO spectrum from 5 to 15 /rm is 
plotted after removal of H 2 O ice and silicate absorption (Boogert et ah 2008). The Gaussian 
at 11.3 /mi (red) corresponds to the PAH detection reported by Bregman et al. (2000). The 
pyrene spectrum (green/blue) was scaled so that the area of band A is equivalent to the 
area of the corresponding YSO feature. All observational and experimental spectra were 
convolved to a matching resolution of ~ 8 cm“^. Vertical lines denote the 5 to 8 /tm region 
where we estimated the PAH contribution to observed absorption. 

The average PAH molecule in the ISM contains 50 carbon atoms (Tielens 2008), rather 
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than the 16 of pyrene. To account for this we assume that the band strengths per bond 
remain constant amongst different PAH species, although computational spectra predict the 
CH stretching band to increase for large PAHs (Bauschlicher et ah 2008), and we report 
contribution percentages in terms of the ratio of CC to CH bonds which characterizes the 
actual mixture of observed interstellar PAHs. For reference, that ratio is = 1.9 for pyrene. 
The simplest PAH, naphthalene (CioHg), has a ratio of = 1.375, and circumcoronene 
(C54H18), the most compact and radially symmetric PAH with ~ 50 carbon atoms, has a 
ratio of = 4. Equations (5) — (7) spell out how to calculate fpAH- In these equations. Ax 
refers to the FWHM of a particular feature. Band strengths from bands B— H in Table 2 
were summed and divided by 19 (the number of CC bonds in pyrene) to obtain a band 
strength of Acc = 4.497 (±0.234) x 10“^® cm per CC bond. The band strength of band A is 
Ach = 6.036 (±0.312) x 10“^® cm per CH bond. 


N{PAH) = 

AcnncH Accncc 

E5_8^„,tAx = (6) 

Ach ncR 

/pah = X 100 (7) 


4.3. Carbon Budget 

Observations of PAH emission from the diffuse ISM indicate that 10 to 20% of the 
total carbon is locked in PAH molecules (Peeters 2011). The column densities of CH bonds 
reported in Table 5 can be used to calculate the fraction of the carbon abundance locked up 
in PAHs in the ices around YSOs. Hydrogen column densities reported by Boogert et al. 
(2013) were converted to carbon column densities using x(C) = 3.9 x 10“'^ (Tielens 2008). 
The number of carbon atoms locked in PAHs were calculated from A(CH) assuming that 
the average PAH molecule contains ~ 22 CH bonds. That value was determined using 
Figure 12 assuming the average PAH contains 50 carbon atoms (Tielens 2008). Table 6 
shows the estimates for the percentages of carbon locked in neutral PAHs around YSOs. 
Sellgren et al. (1995) estimated this value for Mon R2 IRS 3 to be about 4% which is half of 
the value reported in this work, but they also used a higher band strength for the CH bond 
of 1.7 X 10“^® cm/bond. With a larger sample of YSOs, Brooke et al. (1999) estimate the 
percentage of carbon in PAHs to range from 8 to 23%. 
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Fig. 12. — Relationship between number of CH bonds and number of carbon atoms for 
PAHs ranging in size from C14H10 through C50H22, as listed in the experimental section 
on www.astrochem.org (Bauschlicher et ah 2010; Boersma et al. 2014). According to this 
relationship, an average-sized PAH containing 50 C atoms (Tielens 2008) should have roughly 
22 CH bonds. 

We chose the CH stretching bands to constrain PAH abundances because they tend to 
be most consistent across PAH species, but Keane et al. (2001) estimated PAH abundances 
toward YSOs using excess absorption at 6.2 pm instead. Their estimates for the percentage 
of total carbon locked in PAHs ranged from 83%, assuming neutral PAHs, to 7%, assuming 
ionized PAHs. This vast difference is due to the fact that neutral PAHs have relatively strong 
3.25 pm features and weak 5 to 10 pm features, while ionized PAHs exhibit the opposite 
behavior. According to our estimates, Mon R2 IRS 3 and S140 IRS 1 have less than 10% 
of their total carbon in neutral PAHs. These estimates are upper limits in terms of neutral 
PAH abundances since they assumed the entire 3.25 pm features were due solely to neutral 
PAHs. If the environments near YSOs follow the same trends as the rest of the galaxy, 
namely that PAHs in these systems contain 10 to 20% of the carbon reservoir, ionized PAHs 
should account for the remaining PAH contribution to the carbon budget. When ionized 
PAHs are accounted for, the portion of absorption in the 5 to 8 pm spectral region due to 
PAHs should be even higher than what is reported in Table 5. The estimate for GL 2136 
that neutral PAHs account for 22% of its total carbon seems high, lending support to the 
conclusion of Brooke et al. (1999) that absorption from crystalline H2O ice is responsible for 


the apparent 3.2 /im “feature” in its spectrum. Therefore, caution must be exercised when 
interpreting the cause of absorption in this region. 


5. Summary and Future Work 

PAHs account for 10 to 20% of the Milky Way’s carbon reservoir, yet their direct 
observation via absorption in the dense interstellar medium remains difficult. Quantification 
of PAHs expected to be embedded in interstellar ices is hindered by the limited knowledge 
of the absolute band strengths, peak wavelengths, and band widths for PAHs frozen in ice 
matrices. Our work addresses this issue via an investigation of pyrene embedded in water 
ice. 

1. We report the hrst laboratory determination of the band strength for the CH stretching 
mode of pyrene in water ice near 3.25 pm. 

2. We report new absolute band strengths from 3 to 15 pm for pyrene embedded in H 2 O 
and D 2 O ice. The band strengths reported here are roughly 50% greater than those 
published by Bouwman et ah (2011). 

3. We combed the dataset of YSO spectra published by Boogert et al. (2008) to look 
for evidence of PAH absorption at 3.25 pm and used the results from our laboratory 
measurements to estimate PAH column densities where applicable. In our sample, 
neutral PAHs account for 2 to 9% of the non-H20 ice absorption from 5 to 8 pm. 
Neutral PAHs account for 5 to 9% of the carbon budget toward YSOs, depending on 
the particular PAH mixture present. 

Further laboratory work is needed to determine widths, peak positions, and band 
strengths for a large variety of PAH species embedded in ices. Figure 11 shows why this 
is critical. It is obvious from a comparison of the observed absorption of Mon R2 IRS 3 
to our laboratory measurements that pyrene alone cannot explain the feature at 3052 cm“^ 
(3.25 pm) attributed to CH stretching. Band A in our laboratory spectra is too narrow 
and slightly offset from the observed 3.25 pm feature. Beyond 10 pm, the CH out of plane 
bending features for pyrene do not align with the PAH detection at 11.3 pm reported by 
Bregman et al. (2000). In addition, the pyrene features in this region are much stronger 
than the reported 11.3 pm feature. This apparent discrepancy supports the conclusion that 
PAHs around Mon R2 IRS 3 must contain a mixture of species. As illustrated in Figure 2 of 
Allamandola et al. (1999), a mixture of neutral PAHs will exhibit a strong, narrow feature at 
3.25 pm and much broader, weaker features beyond 10 pm. This is because the position of the 
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Table 5. Percent Contribution to the 5 to 8 /rm Region (tpAn) After H 2 O Subtraction 


Object 

iV(CH) (10^® cm-2) 

iV(CC) (10 

cm ^) 

fpAH (%) 



^ = 1.375 

ncn 

^cc _ 4 
non 

^ = 1.375 

_ 4 

ncn 

Mon R2 IRS 3 

3.76 ± 0.41 

5.17 ± 0.56 

15.03 ± 1.62 

2.9 

8.5 

S140 IRS 1 

2.05 ± 0.23 

2.81 ± 0.32 

8.18 ± 0.92 

2.0 

6.0 

GL 2136*" 

13.34 ± 2.30 

18.34 ± 3.16 

53.36 ± 9.18 

6.9 

20.4 


Note. — ^^Brooke et al. (1999) report that the inflection point of a crystalline ice absorption 
profile coincides with the apparent “feature” in GL 2136, so the calculated abundances reported 
here for this system may be inaccurate. 


Table 6. Carbon Budget for Neutral PAHs Around YSOs 


Object 

Ng 

(10^^ cm“^) 

Nq, total 
(10^® cm“^) 

iV(GH) 
(10^® cm“^) 

Nq, in PAHs 
(10^® cm“^) 

% Garbon 
in Neutral PAHs 

Mon R2 IRS 3 

26.19 

102.14 

3.76 

8.55 

8.37 

S140 IRS 1 

22.41 

87.40 

2.05 

4.66 

5.33 

GL 2136*^ 

35.20 

137.67 

13.34 

30.32 

22.02 


Note. — ’^Reported by Boogert et al. (2013). ^Brooke et al. (1999) report that the inflection 
point of a crystalline ice absorption profile coincides with the apparent “feature” in GL 2136, so 
the calculated abundances reported here for this system may be inaccurate. 
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CH stretching mode is quite stable across PAH species whereas the CH out of plane bending 
modes, though stronger for any individual species, are much more variable in position. The 
result is an apparent broadening and weakening of the composite absorption feature beyond 
10 /im relative to the 3.25 /im feature in the spectrum of a PAH mixture. 

Future work should also investigate the effects of other relevant ice matrices (e.g. CO 2 ) 
on the PAH absorption profiles. Eventually, the infrared band strengths for a variety of 
ionized PAHs need to be measured to compliment work with neutral species. Finally, 
the greatest difficulty for detecting PAH absorption features in YSO spectra is the lack 
of sufficient signal-to-noise due to atmospheric effects (in particular telluric CH4 lines) in 
the ground-based data. Observational campaigns during the upcoming SOFIA and JWST 
missions should effectively solve this problem. 
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